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A new approach to separation of a variety of microscopic and mesoscopic objects in dilute solution 
is presented. The approach takes advantage of unique properties of a specially designed separation 
device (sieve), which can be readily built using already developed microlithographic techniques 
Due to the broken reflection symmetry in its design, the direction of motion of an object in the 
sieve varies as a function of its self-diffusion constant, causing separation transverse to its direction 
of motion. This gives the device some significant and unique advantages over existing fractionation 
methods based on centrifugation and electrophoresis. 

PACS numbers:82.45.+z, 87.f5.-v, 07.f0.-h, 36.20.Ey 



Separation of macromolecules such as proteins and 
DNA, as well as mesoscopic objects such as cells and latex 
spheres according to size has many important technolog- 
ical applications, and an immense effort has gone into 
achieving efficient, well controlled and high resolution 
separation techniques |l| [| . Two of the main avenues 
that have been employed extensively are electrophoresis 
Pj and centrifugation f|. These two methods are some- 
what complementary since the former typically separates 
by charge density and the latter by mass density. In par- 
ticular, virtually all electrophoretic techniques rely on an 
electrophoretic mobility /i e that changes as a function of 
molecular weight M w , or some other characteristic for 
which separation is desired, since the separation occurs 
along the same direction as the average motion. An ini- 
tially polydisperse band separates into many bands con- 
taining objects of different sizes as they travel at different 
velocities v(M w ) = ^ e (A/ w )E along the direction of the 
applied field E. A major obstacle to the electrophoretic 
separation of large polyelectrolytes with constant charge 
density such as nucleic acids is the independence of their 
electrophoretic mobility to their molecular weight in so- 
lution. To achieve separation, the polyelectrolytes are 
typically placed in a gel medium, where steric interac- 
tions generate a size-dependent /i e - Despite significant 
progress in refined gel electrophoresis techniques, large 
objects such as cells or subcellular structures are impos- 
sible to separate due to the limited range of achievable 
pore sizes, and issues such as sample loading and recovery 
are especially problematic for fragile specimens j^]. 

In this paper, a completely new approach to separa- 
tion, which embodies the advantages of both free flow 
electrophoresis and gel electrophoresis, and is made 
possible by microlithographic techniques recently intro- 
duced by Volkmuth and Austin 0, is presented. The 
general idea is to design a particular electrophoresis 
chamber (sieve) such that the electrophoretic mobility 
tensor is non-diagonal, i. e., the objects do not move 
along the direction of applied electric field E. Further- 
more, the direction of motion varies as a function of size, 



causing objects of different sizes to move along different 
directions. As a result, this sieve causes lateral separation 
as in free flow electrophoresis |I) without the complica- 
tions involved in creating a uniform transverse laminar 
flow field. The sieve is reusable, and the approach en- 
ables continuous operation since the paths are separated 
spatially rather than temporally, making retrieval of sep- 
arated products extremely easy. These properties po- 
tentially provide a significant advantage over traditional 
methods for separation of large quantities. 

In the remainder of this Letter, a specific realization 
of such a sieve is presented. The dynamics of a poly- 
electrolyte, e. g., a DNA molecule, subject to a uniform 
electric field in the sieve, and the mobility and resulting 
separation characteristics are derived. These results are 
then verified by numerical simulation, followed by a de- 
tailed comparison of the method with existing separation 
techniques. 
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FIG. I. The geometry of the sieve. Molecules enter the top 
of a cell from a narrow opening and diffuse away from the left 
wall as they are pulled down by the electric field. Smaller 
molecules diffuse farther and are therefore more likely to end 
up to the right of the branching point, located at a distance 
yb from the entrance and Xb from the left wall. Therefore, the 
branching probability pb, and subsequently the macroscopic 
mobility of a molecule, depends on its size. 
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The geometry of the sieve is shown in Fig. [|. It con- 
sists of a rectangular array of "cells" of periodicity x X yo 
that have a narrow entrance at the top and two exits 
at the bottom, which connect to the next row of cells. 
DNA fragments, characterized by a persistence length 
£ p , contour length L — N£ p , and diameter d, enter from 
the top and move down the cell subject to an electric 
field E = Ee y . The radius of gyration R g and the self- 
diffusion constant D G of a DNA fragment of N persis- 
tence lengths are given by 
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D G ~ D N~ a , 
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where v = 1/2 is the swelling exponent, Do is the diffu- 
sion constant for a single persistence length, and a = 1 
for the Rouse model ||] (free draining conditions) and 
a = v = 1/2 for the Zimm model H, where hydrody- 
namic interactions are taken into account 0. The clcc- 
trophoretic mobility of the fragments is independent of 
size and given by 



A 

Me = J, 



(3) 



where A and £ are the effective charge and friction coef- 
ficient of a persistence length of the DNA, respectively. 
(The Reynolds number is very small in all cases of in- 
terest, and inertial effects can be ignored.) Note that 
hydrodynamic interactions induced by Brownian Motion 
are not screened by counter-ions, as is the case of elec- 
trophoretic velocity fields |p|^0||. Consider a fragment 
that enters a cell at time to and diffuses away from the 
left wall as it drifts down the cell. Ignoring its internal 
modes and characterizing its dynamics simply by its elec- 
trophoretic mobility fi e and self-diffusion constant D G , 
the probability of finding the DNA a distance x away 
from the left wall at time t when it has drifted a distance 
y{t) = /J, e E(t — to) from the top of the cell is given by 
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which can be obtained from the solution to the diffusion 
equation with reflective (Dirichlet) boundary conditions 
on the left wall, and ignoring the effect of the right wall. 
A branching point is located at a distance Xb from the 
left wall and y& from the top of the cell. Fragments that 
have diffused farther from the left wall than the branch- 
ing point are funneled to the entrance of the cell diagonal 
to the original one, whereas those that are closer to the 
wall end up at the entrance to the cell immediately be- 
low. The probability of branching can be calculated from 
Eq.(|) as 
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For a DNA fragment of N persistence lengths, 

PB^-ie-W", (6) 



where 



N 



(7) 



is the characteristic separation size of the sieve, and A 
is a constant fitting parameter of 0(1). A remarkable 
observation is that No can be tuned simply by changing 
the applied electric field, increasing the dynamic range of 
separation dramatically. 

A more accurate analytical estimate of pb (A) requires 
taking into account internal relaxation of the polymer, 
effective wall potentials at distances less than the ra- 
dius of gyration of the polymer, corrections due to dif- 
fusion along the field direction, among other effects. For 
example, the right wall can no longer be ignored for 

x ° ) , and the finite size of the fragment 

becomes significant when R g > Xb, or equivalently N > 
(xb/ip) 1 ^ ■ Relaxation effects will be important when 
cell traversal time yb/(^ e E) is less than the principal re- 

laxation time Ri/D G , i. e., for N > ( — - 

9 \ D y b J 

Computation of these effects is beyond the scope of this 
Letter. However, from a practical standpoint, characteri- 
zation for a given geometry can be more readily achieved 
by numerical simulation or experiment for most appli- 
cations, since pb(N) completely characterizes the sep- 
arational properties of the sieve and the periodic struc- 
ture greatly simplifies the numerical determination of this 
quantity by simulating a single cell. 

At this point, it is useful to point out that any sieve 
design that breaks the left-right symmetry can in prin- 
ciple be used for purposes of separation. The underly- 
ing idea is reminiscent of rachet potentials [pl]Jl^ 1 , which 
have also been recently proposed as particle separators 
HQ. Both rachets and the sieve exploit differences 
in diffusion constants and steric constraints, and if the 
y axis is identified with time, the sieve can actually be 
thought of as the time history of a rachet potential. How- 
ever, unlike rachets, the sieve does not require time de- 
pendent potentials and is able to operate in continuous 
mode. 

Numerical integration of the equations of motion using 
the Rouse Model for polymers ||J^] (a = 1) and reflect- 
ing boundary conditions at the walls is in good agreement 
with analytical results. Figure |^ shows how the proba- 
bility density for the COM position of a polymer with 
N = 24 evolves as it moves down the cell. After the ini- 
tial internal relaxation time, the probability can be fitted 
to the form of Eq.(||); the inset shows the expected lin- 
ear increase in the mean square distance (x 2 ) from the 
left wall as a function of distance y from the top of the 
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cell. Figure || shows the branching probability p B (N) as 
a function of polymer size, and has the expected expo- 
nential behavior, even though the sizes of the polymers 
become quite large compared to the cell feature size. A 
correlation analysis of the time series of branching events 
confirms that they are statistically independent from one 
cell to the next. 



0.8 




10 20 

distance from wall x 

FIG. 2. Evolution of the density profile of chains with 
N = 24 persistence lengths during the traversal of a cell. 
After an initial relaxation period, the profile approximately 
takes the form given by Eq. (|B|), and the mean square dis- 
tance from the left wall increases linearly as a function of the 
distance from the entrance, as demonstrated in the inset. In 
this run, the primitive cell size of the rectangular lattice is 
xo = 40 x j/o = 80; yt = 60 and xt = 4.8 (all in arbitrary 
units.). 
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FIG. 3. Branching probability as a function of the number 
of segments N. The line is a fit to the exponential form (^) 
with v — 1 and TVo ~ 49. 

Although this method of separation sounds very 
promising in principle, it is important to assess per- 
formance parameters and feasibility before a decision 
can be made about its practicality and whether it 
can compete against established techniques for cer- 
tain tasks. One of the most important issues is the 
resolution that can be achieved Jlal. A monodis- 



perse packet of polymers with size N will spread lat- 
erally as it moves through the sieve, and after passing 
through M rows, the density profile of te band will ex- 
hibit a Bernoulli distribution whose peak is located at 
Xp(M) — xqMpb(N) and whose variance is <J 2 {N) — 
xlM PB (N)[l - p B (N) @. Hence, the full width at 
half maximum (FWHM) of the corresponding band will 
be FWHM(A) ~ 2x^2\n(2)p B {N)[l -p B (N)}M. On 
the other hand, the peak separation between polymers of 
sizes N and N+SN in a polydisperse sample will increase 
as SXp(N) = Mx (SN)[dp B /dN}. Thus, resolution can 
be improved indefinitely by passing the polymers through 
more rows of cells. Optimal resolution is achieved when 
N ~ No, for which M ~ (5N/N) 2 rows are needed in or- 
der to resolve these two peaks. For a cell size of about 5 
microns, up to M — 10 4 rows should be feasible, enabling 
single-segment resolution up to N ~ 100 and 1% resolu- 
tion beyond that. Note, however, that resolution can be 
further enhanced by gradient methods that are frequently 
implemented in gel electrophoresis [[ItJ , in this case by a 
spatially varying electric field or cell size. 

Another major concern is the throughput, which is af- 
fected by various factors, including cross-sectional area, 
concentration and velocity of the polymers, and ease of 
specimen loading and extraction. Original electrolitho- 
graphic designs proposed and tested by Austin and 
coworkers were extremely shallow |[],[l8), with depths 
comparable to £ p , in order to maximize hooking and to 
enable individual visualization of polymers. The sieve de- 
sign presented here, on the other hand, will benefit from 
increased depth to at least R g , and the cells should be 
designed as deep as possible to increase throughput. A 
stacked configuration might be considered if the mechan- 
ical stability of posts becomes a concern. A major bottle- 
neck is the entrance to the sieve, since all polymers should 
start from the same point, rather than a band, to achieve 
separation. This will give rise to significantly increased 
concentrations near the entrance to the sieve. Although 
a dilute solution is assumed in the calculation presented 
here, separation is not limited to the dilute regime and 
occurs in semidilute solutions as well. Since the mobility 
of polymers is significantly higher in the absence of a gel 
medium, much higher concentrations can be tolerated, 
therefore this bottleneck may not be as problematic as 
it appears. Furthermore, the easy and quick extraction 
of the separated specimens enables continuous operation, 
in which polymers are constantly added at the entrance 
and extracted from exit channels placed at the bottom 
of the sieve. Although the sieves might be expensive to 
produce, they are tunable and reusable, significantly low- 
ering their effective cost. Thus, this technique may have 
major advantages over traditional methods for separation 
of large quantities. 

The same technique can be applied to separation 
through centrifugation pi as well. Individual sieves 
shaped as pie slices, with their entrances at the apex, 
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can be arranged as stacked pies. The polydisperse solu- 
tion can then be fed through a tube along the rotation 
axis and separated objects can be collected in contain- 
ers around the circumference that rotate with the sieves. 
Virtually any separation characteristic can be achieved 
by modulating the cell size as a function of distance to 
the rotation axis. 

The particular design presented here is only one possi- 
bility and is by no means an optimal geometry, although 
its performance is expected to be quite satisfactory. On 
the other hand, the scalable structure of the sieve makes 
it possible to fully characterize device performance from 
a detailed modeling of dynamics within a single cell, sig- 
nificantly simplifying a numerical design effort. Further- 
more, technological requirements for an experimental re- 
alization is well within today's capabilities, as has been 
demonstrated by Volkmuth et al. Therefore, exper- 
imental verification of the soundness of the general ap- 
proach and a feasibility study of a prototype device is 
within reach in the very near future. 
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